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Pref ace 

This  thesis  is  in  response  to  the  continued  interest  in 
the  PRIZ  device  at  AFIT.  Some  experimental  work  has  been 
done  in  the  United  States,  primarily  at  AFIT,  but  the 
majority  of  the  theoretical  work  to  date  has  been  in  the 
Soviet  Union.  A  review  of  the  theoretical  work  done  in  the 
Soviet  Union  is  presented  here,  and  further  development  of 
the  theory  is  presented.  A  numerical  model  of  the  PRIZ 
device  is  presented  and  compared  with  theory  and  experiment. 

My  sincerest  thanks  go  to  Or.  Bailey  for  his  guiding 
hand  and  patience  throughout  this  project.  Perhaps  he  will 
be  rewarded  with  a  competent  researcher.  Dr.  Luke  and  Major 
Lupo  deserve  special  mention  for  providing  their  inputs  into 
this  project  as  well. 

Special  thanks  is  due  Major  Lupo,  whose  graphical 
subroutines  added  immeasurably  to  this  thesis. 

Finally,  I  would  like  to  thank  my  family,  who, 
although  their  support  was  of  a  long  distance  nature,  always 
deserve  at  least  some  of  the  credit  for  anything  I  do. 

Xav^Ld.  A .  Sw&Jt-Lng- 
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Abstract 

A  numerical  computer  model  of  the  PRIZ  electro-<optic 
spatial  light  modulator  is  presented.  Quided  by  the  earlier 
work  of  Astratov,  this  model  involves  differencing  the 
longitudinal  spatial  derivative  and  using  a  numerical 
integration  technique  to  solve  the  time  derivative.  The 
model  is  validated  by  comparing  the  numerical  results  with 
the  analytical  theory , developed  by  Bryskin.  Bryskin  has 
spent  several  years  developing  a  theoretical  analysis  of  the 
PRIZ  device  (ref.  10  and  ref.  9).  The  primary  elements  were 
reviewed  and  in  selected  instances ,  extended.  Differences 
are  noted  and  possible  explanations  are  offered. 

'Comparisons  are  also,  made  with  experimental  results. 

Finally,  the  model  is  extended  to  solve  for  the  integrated 
transverse  electric  field  in  the  case  of  the  simple  spot. 

The  integrated  field  is  found  to  have  nearly  a  1/r 
dependence.  The  solution  for  a  simple  spot  is  used  to 
explain  the  "edge-enhancement'’’  effect  of  a  simple  spot 
reported  by  Shields.  The  damage  to  the  PRIZ  is  then  briefly 
discussed. 
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Chapter  I 


PRIZ  is  a  Russian  acronym  for  a  spatial  light 
modulator.  The  PRIZ  device  is  a  simple  modulator  capable  of 
optical  image  storage  and  even  dynamic  imaging  with  edge 
enhancement.  Potential  applications  in  the  areas  of  optical 
computing,  optical  signal  processing,  and  image  pattern 
recognition  (12:1)  have  stimulated  significant  research  and 
development  efforts  on  these  devices.  Petrov  and  Bryskin 
(ref.  10  and  ref.  16)  have  led  a  concerted  and  sustained 
Soviet  developmental  effort  which  has  provided  experimental 
data  on  the  PRIZ.  They  have  formed  a  theoretical  basis  for 
understanding  the  principles  of  operation  and  have 
characterized  the  modes  of  operation.  Experimental  results 
consist  chiefly  of  longitudinal  electric  field  measurements 
in  various  modes.  The  theoretical  development  has  provided 
analytic  expressions  for  certain  limiting  cases,  and 
highlighted  two  operational  regimes  of  the  device. 

Recently,  Shields  and  Luke  (ref.  17)  constructed  and 
tested  a  PRIZ.  Their  studies  (the  first  outside  of  the 
Soviet  Union),  confirmed  those  capabilities  reported  by  the 
Soviets,  and  reproduced  many  of  their  results.  In  a 
subsequent  study,  Nilius  and  Luke  (ref.  12)  reported  on  the 
memory  capability  of  the  device.  Anderson  and  Luke  (2:3) 
later  noted  that  under  typical  "write  beam"  conditions,  the 
crystal  can  experience  significant  damage.  This  damage  can 
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lead  to  modification  and  degradation  in  the  PRIZ  operating 
characteristics;  thus  severely  restricting  development. 

Problem  Statement 

A  generalized  theoretical  analysis  was  initiated  to 
parameterize  the  operating  regime  of  the  PRIZ  devices. 
Previous  analyses  provided  guideposts  in  the  validation 
process;  but  were  of  limited  or  uncertain  general  utility. 
They  often  employed  unrealistic  or  restrictive  assumptions 
which  were  inconsistent  with  the  conditions 
encountered  for  realistic  device  applications.  The 
generalized  approach,  however,  was  baaed  on  a  direct 
numerical  integration  of  the  rate  equations  governing  FRIZ 
operation.  The  kinetic  processes,  either  in  number,  spatial 
dependence  or  type,  were  not  restrictive  to  the  method 
employed.  Employing  this  approach,  the  temporal  and  spatial 
distribution  of  charges  and  the  internal  fields  have  been 
characterized  under  a  broad  range  of  excitation  conditions. 
Comparisons  have  been  made  with  existing  analytic  solutions, 
and  scaling  laws  have  been  identified. 

PRIZ  Fabrication  and  Operation 

The  PRIZ  is  most  often  built  of  the  photo-refractive 
crystal,  Bi^SiOgQ  (BSO)  (17:1).  Some  experiments  relating 
to  the  PRIZ  have  been  performed  with  the  crystal,  Bij2f'eO20 
(3:1585).  These  crystals  form  the  basis  for  two  types  of 


2 


devices:  the  'conducting’  PRIZ  and  the  ’standard’  PR  I Z . 

These  are  pictured  in  Figure  1.1.  The  ’conducting’  PRIZ 
consists  of  s  (111)  or  (110)  cut  BSO  crystal  to  which 
trsnsparent  electrodes  are  applied.  The  ’standard’  PRTZ 
consists  of  the  sene  BSO  crystal;  however  in  this  design  the 
crystal  is  sheathed  in  dielectric  aaterial  before  the 
application  of  the  transparent  electrodes. 

The  external  electric  field  is  oriented  parallel  to  the 
incoming  write  beaa.  The  write  bean  typically  deposits 
anywhere  froa  1  to  2S0  w J  on  the  device.  The  wavelengths 
used  are  in  .he  visible  range,  so  the  device  can  be  viewed 
without  sophisticated  optical  equipment.  The  write  b**a 
photo-exc i tea  electrons  froa  donor  sites  into  shallow  traps, 
between  0.6  eV  and  0.3  eV  (10:878).  Froa  the  shallow  traps, 
electrons  can  be  theraally  excited  into  the  conduction  band . 
Write  beaaa  coaposed  of  aore  energetic  photons  an  e\  it* 
electrons  directly  into  the  conduction  band.  The  ele<  t  r  r.s 
in  the  conduction  band  drift  due  to  the  external  appi i ed 
field  and  are  trapped  in  deeper  traps,  at  I . 1  *i  and  .  .5  - . 
(15:3686).  The  trapped  electrons,  the  now  positive), 
charged  donor  sites,  and  the  still  aobile  electrons  rente  » 
non-unifora  charge  distribution  m  the  crystal.  T his  -barge 
distribution  acts  to  set  up  a  new  internal  field  with 
transverse  coaponenti. 

Using  a  polarizer  before  the  tlluainated  face  and  an 
analyzer  after  the  dark  face,  Figure  1.2.  a  resd  bea«  an  be 
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gure  1.2.  Complete  PRIZ  with  polarizers  for  readout 
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sent  through  the  device  to  interrogate  the  optical  changes 
induced  by  the  charge  distribution  in  the  crystal.  The 
charge  distribution  is  a  map  of  the  write  beam  which  in  a 
full  application  of  the  device  would  result  from  some  kind 
of  image.  The  read  beam  is  sent  through  the  device, 
parallel  to  the  write  beam.  It  is  of  a  longer  wavelength 
than  the  write  beam.  By  havihg  the  read  beam  at  a  longer 
,  wavelength,  the  absorption  coefficient  is  generally  lower 

(17:12),  as  is  the  quantum  efficiency  (17:12).  This  reduces 
the  number  of  electrons  photo-excited  by  the  read  beam.  The 
internal  charge  distribution  of  the  crystal  and  the 
resulting  electric  field  serve  to  modulate  the  index  of 
refraction  via  the  linear  electrooptic  effect  (14:646). 

Since  the  internal  electric  field  of  the  crystal  is 
spatially  varying,  the  local  polarization  of  the  beam  is 
also  spatially  varying.  The  analyzer  then  provides  a 
corresponding  intensity  variation  (14:646). 

Now  that  we  have  established  the  origin,  nature,  and 
importance  of  the  internal  electric  field  we  can  follow  the 
development  of  Tanguay  (14:646)  to  see  how  the  phase 
difference  of  the  device  with  a  (111)  crystal  orientation 
can  be  calculated.  The  phase  difference,  r,  is  expressed  as 
(16:817). 


(1.1) 
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where  X  is  the  wavelength  of  the  incident  light,  n’  is  the 
unperturbed  index  of  refraction  and  dQ  is  the  longitudinal 
dimension  of  the  crystal.  ,  the  integrated  transverse 

component,  is  given  by 


E 


t 


Etrans<x^'z) 


dz 


(1.2) 


where  Etrang  is  the  electric  field  component  perpendicular 
to  the  z  axis. 

The  device  can  now  be  treated  with  standard 
electrooptic  formalism,  provided  the  phase  difference  is 
known.  The  quantity  presenting  the  most  difficulty  in 
determining  the  phase  difference  is  the  longitudinally 
integrated  transverse  electric  field.  The  balance  of  this 
paper  will  be  devoted  to  finding  a  solution  in  some  cases 
for  this  integral.  The  longitudinal  field  and  the 
longitudinal  charge  density  can  be  used  to  determine  the 
transverse  field. 


Sequence  of  Presentation 

Chapter  II  introduces  the  continuity  and  charge 
transport  equations  used  to  model  the  PRIZ  device.  The 
boundary  conditions  affecting  the  mathematical  solution  of 
the  basic  equations  are  discussed,  and  related  to  the 
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physical  structure  of  the  device.  Existing  analytic 
solutions  are  then  presented. 

The  method  of  solution  used  in  the  computer  model  is 
presented  in  Chapter  III.  The  numerical  techniques  used  in 
the  solution  are  reviewed,  and  the  results  of  the  model  are 
presented.  Where  applicable,  the  solutions  generated  by  the 
model  are  compared  to  predictions  taken  from  the  analytic 
solutions  developed  in  Chapter  II.  The  model  is  also 
extended  into  realms  into  which  the  analytic  solutions  are 
inapplicable . 

In  order  to  qualitatively  evaluate  PRIZ  operation, 
including  the  transverse  field,  a  simple  example  of  an 
illuminated  spot  is  treated  in  Chapter  IV.  The  potential 
distribution  through  out  the  crystal  is  calculated  and  used 
to  solve  for  the  transverse  fields,  which  are  then 
integrated  to  give  a  solution  for  Et*.  Previous 
experimental  observations  are  explained  using  this  model. 

Chapter  V  summarizes  the  results  and  includes 
recommendations  for  further  theoretical  and  experimental 
efforts . 
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Chapter  II 


In  order  to  investigate  the  temporal  evolution  of  the 
fields  within  the  PRIZ,  the  excitation  process  and 
subsequent  charge  transport  and  trapping  are  analyzed.  A 
qualitative  discussion  of  the  mechanisms  involved  is 
presented,  followed  by  an  analysis  of  the  basic  equations  of 
generation  and  transport.  The  relation  of  the  type  of  PRIZ 
to  the  specifications  of  boundary  conditions  on  the 
solutions  is  then  examined.  Finally,  analytic  solutions  in 
several  limiting  cases  are  presented. 

Qualitative  Discussion 

Before  examining  the  basic  equations  governing  PRIZ 
operation  and  identifying  the  operating  regimes,  it  is 
helpful  to  discuss  the  PRIZ  from  a  qualitative  standpoint. 
The  PRIZ  operates  by  using  an  incident  light  beam  to 
photo-excite  electrons  into  the  conduction  band.  Figure  2.1 
is  a  diagram  of  the  band  structure.  The  primary  source  of 
electrons  is  the  Si  vacancy  at  2.6  eV.  The  luminescence 
centers  at  1.3  and  2.25  eV  can  also  contribute  electrons  to 
the  conduction  band,  if  they  are  occupied.  The  excitation 
rate  out  of  these  levels  is  proportional  to  the  power  of  the 
incident  light.  The  absorption  of  the  incident  light  in  the 
bulk  of  the  sample  results  in  the  excitation  rate  displaying 
an  exponential  decay  consistent  with  Beer’s  law.  Electrons 
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Figure  2.1.  Energy  Band  Diagram  of  Bismuth  Silicon 
Dioxide  at  Room  Temperature  (15:3684) 


excited  out  of  these  levels  leave  behind  a  positively 
charged  donor  site. 

Once  in  the  conduction  band,  the  electrons  can  then 
decay  back  into  any  one  of  the  levels  that  has  a  hole 
available  for  recombination.  The  probability  of  the 
electron  falling  into  the  level  is  determined  by  the  density 
of  available  sites  to  decay  in,  and  the  depth  of  the  level 
in  the  band  gap.  Electrons  which  decay  into  the  shallow 
traps  can  be  reexcited  with  relative  ease  into  the 
conduction  band.  The  level  at  1.3  eV  is  called  a  deep  trap, 
due  to  the  fact  that  electrons  cannot  be  thermally  excited 
out  of  this  level,  and  the  write  light  under  consideration 
(400  to  500  nm. )  does  not  have  the  appropriate  energy  to 
excite  these  levels  to  the  conduction  band.  Thus  conduction 
band  electrons  recombining  to  this  level  result  in  trapped, 
negatively  charged  centers. 

The  externally  applied  electric  field  causes  transport 
of  the  electrons  injected  into  the  conduction  band.  For  the 
standard  bias  conditions,  the  electrons  drift  away  from  the 
illuminated  electrode  and  the  holes  drift  toward  it.  Due  to 
the  relative  immobility  of  the  holes  for  the  times  under 
consideration,  they  will  be  treated  as  fixed  throughout  this 
study.  Since  the  holes  are  immobile,  and  £he  electrons  are 

moving  away  from  the  illuminated  electrode,  a  region  of 
positive  space  charge  is  formed  just  beyond  the  surface  of 
the  illuminated  electrode.  The  shallow  traps  at  0.3  eV  and 
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0.6  eV  (10:878)  can  trap  the  electrons.  Subsequent  thermal 
excitation  leads  to  a  reintroduction  of  the  electrons  into 
the  conduction  band.  The  net  effect  is  an  apparent  reduction 
in  the  mobility  of  the  electrons  in  the  conduction  band.  In 
contrast,  since  the  electrons  are  trapped  relatively 
permanently  by  the  1.3  eV  level,  a  region  of  negative  charge 
can  form  beyond  the  positive  charge  region.  This  process  of 
generation,  transport  and  subsequent  trapping  leads  to 
charge  separation.  Between  these  two  regions  exists  an 
interface  where  the  charge  density  is  neutral.  The  location 
of  the  interface  is  dependent  on  the  ratio  of  the  electron 
drift  length  and  the  absorption  length. 

When  the  total  number  of  electrons  excited  into  the 
conduction  band  is  small,  the  charge  densities  created  in 
the  crystal  do  not  generate  a  large  enough  self  field  to 
materially  effect  the  applied  field.  This  is  referred  to  as 
the  linear  regime.  The  total  longitudinal  electric  field  in 
this  regime  can  be  approximated  by  the  applied  electric 
field.  The  linear  regime  is  characterized  by  a  stationary 
interface . 

As  the  total  number  of  electrons  excited  into  the 
conduction  band  increases,  the  self  field  generated  by  the 
charge  separation  becomes  significant  cpmpared  to  the 
applied  field.  The  total  charge  shielding  the  applied  field 
is  greatest  at  the  interface,  thus  the  self  field  is 
largest  at  that  point.  Since  the  self  field  is  oriented  in 
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the  opposite  direction  of  the  applied  field,  the  magnitude 
of  the  total  field  is  a  minimum  at  the  interface  point. 

This  is  referred  to  as  the  "bottleneck".  The  velocity  of 
the  electrons  in  the  conduction  band  is  the  product  of  their 
mobility  and  the  local  value  of  the  electric  field.  Since 
at  the  bottleneck  the  field  obtains  its  minimum  value,  the 
electron  velocity  is  also  at  its  minimum.  This  velocity 
reduction  results  in  an  enhanced  trapping  in  the  vie  inity  of 
the  bottleneck,  and  a  resultant  migration  of  the  interface 
point  toward  the  illuminated  electrode.  Thus  in  the 
nonlinear  regime,  the  interface  moves  toward  the  illuminated 
face  of  the  crystal. 


Basic  Equations 

The  continuity  equation  yields  the  temporal  variations 
of  the  local  charge  density  within  the  crystal  (5:82).  The 
relevant  equations  for  the  densities  of  interest  are: 


d 

e — (  n  +  n  -  n 

at 


v  j 


(2.1) 


- -  =  gF ( B) 

dt 


(2.2) 


3n 
d  t 


n(nt  -  n_) 
•»* 


n_ 

Tth 


(2.3) 
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The  convective  transport  is  described  by  the  current  density 
equation . 

j  =  e«(  nHQ  -  nvv  +  (kT/e)Vn  )  (2.4) 

where 

e  =  the  modulus  of  the  electron  charge 

gF(B)  -  rate  of  electron  excitation  into  the  conduction 
band 

R  =  the  position  vector 

n+  s  density  of  positively  charged  donor  sites 
n_  s  density  of  trapped  electrons 
n  s  density  of  electrons  in  the  conduction  band 
nt  s  density  of  traps 
j  s  current  density 

s  recombination  time  at  deep  traps 
=  recombination  time  at  donor  sites 
Tth  ~  thermal  ionization  -time  of  shallow  traps 
v>  -  potential  due  to  internal  charge  distribution 
BQ  -  externally  applied  electric  field 
h  -  electron  mobility 

In  order  to  simplify  this  system  of  equations,  we  will 
initially  eliminate  all  nonlinear  terms.  Further,  we  will 
assume  that  the  rate  of  excitation  is  transversely  uniform, 
thus  allowing  gF(R)  to  be  expressed  as  simply  g  exp(-az), 
with  a  being  the  absorption  coefficient.  The  thermal 
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diffusion  of  charge  in  the  current  expression,  equation  2.4, 
will  be  ignored  as  this  is  small  compared  to  the  motion 
induced  by  the  electric  fields.  Thermal  excitation  out  of 
traps  will  be  ignored.  The  thermal  excitation  is  small 
compared  to  the  photoexcitation.  The  total  electric  field 
will  be  expressed  by 


B  =  B0  -  V<*  (2.5) 


Within  these  approximations,  the  system  of  equations  is: 


d 

— (n  +  n  -  n.)  =  ^v(nB)  (2.6) 

at  ~  + 


dn  + 

-  =  g  exp(-az)  (2.7) 

dt 


an  n 


(2.8) 


where  Tt,  the  lifetime  of  an  electron  in  the  conduction 
band,  is  given  by  0f/nt.  Equation  2.6  relates  the  time  rate 
of  change  of  the  total  charge  density  of  the  sample  to  the 
divergence  of  the  current.  Equation  2.7  states  that  the 
time  rate  of  change  of  the  positively  charged  donor  site 
density  is  simply  the  excitation  rate,  g(z).  While  in 
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equation  2.8,  the  time  rate  of  change  of  the  negatively 
charged  trap  site  density  is  governed  by  the  decay  rate  of 
electrons  from  the  conduction  band. 

We  are  initially  interested  in  only  the  longitudinal 
electric  field  in  equation  2.6  because  we  must  have  the 
complete  solution  for  the  longitudinal  field  and  space 
charge  density  before  examining  the  transverse  fields. 
Substituting  the  equations  2.7  and  2.8  into  equation  2.6 
yields  the  time  rate  of  change  of  the  conduction  band 
density. 

dn  n  d 

. -  =  g(z)  -  —  +  u — (  En  )  (2.9) 

at  rt  dz 

where  E  is  the  z-component  of  the  vector  B  and  z  is  the 
direction  parallel  to  the  light  beam  propagation.  With 
equations  2.7-9  we  have  a  complete  description  of  our 
simplified  system. 

Boundary  Conditions 

Two  different  PRIZ  devices  have  been  constructed  or 
described  previously  in  the  literature  (12:1).  The  first  is 
the  "standard"  PRIZ,  Figure  1.1b.  Here  the  crystal  is 
separated  from  each  electrode  by  a  dielectric  layer.  Thus 
the  crystal  remains  charge  neutral  since  no  charge  can 
traverse  the  dielectric  barriers.  The  boundary  conditions 
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on  the  currents  at  either  face  of  the  crysta.  nr- 


j  (  0  ,  t )  =  0  J  .  ) 

j ( dQ , t )  s  0  _  ’  . 

where  dQ  is  the  length  of  the  crystal  in  the  z  direct;  n 

The  second  PRIZ  device  is  the  "conducting''  PR  I 2 ,  Figure 
1.1a.  The  insulating  dielectric  layer  is  omitted  from  th;i 
device.  Without  the  insulating  dielectric,  the  boundary 
conditions  are  described  by 

j(0,t)  =  injection  current  ^.1. 

j(dQ,t)  =  E(d0,t)  n(dQlt)j  t2.ll 

where  the  injection  current  will  be  dependent  on  the 
electric  field  at  the  illuminated  face  of  the  crystal  versus 
the  barrier  potential.  In  moat  cases  the  injection  current 
represented  by  equation  2.12  is  very  small  compared  to  the 
current  represented  by  equation  2.13,  so  j(0,t)  -  0.  The 
bulk  of  this  paper  will  deal  with  the  "conlucting'  vers;  n 
of  the  PRIZ. 
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4*  1 1  tcusaed  previously  m  the  qualitative  discussion, 
there  are  two  operating  regiaea  of  the  PR  1 2  which  depend  on 
the  total  charge  generated  by  the  write  beaa.  The 
excitation  rate,  g(z),  can  be  related  to  the  intensity,  I  of 
the  write  beaa  by  ( 10:879) 

Illlaf 

gin  =  - a  (2.14) 

hr 


where  f }  i j  the  exper i aenta 1 1 y  obtained  quantua  efficiency 
at  any  given  wavelength  and  hr  is  the  photon  energy.  A 
eharac t er i at i c  tiae,  t  ,  divides  the  linear  and  nonlinear 
regiaea  >f  operation,  where  (9:1349): 


t 


o 


n  d_ 

—a  a 


g(0)afBol rt 


(2.15) 


with 


«  ’ 

nQ  -  -  i  2  .  16) 

•-rt 

where  ♦  ’  is  the  total  dielectric  constant  of  f  he  aaterial. 
Physically,  the  tue  t}  is  the  length  of  t  i  ae  that  aust 
elapse  before  the  self  field  is  the  saae  order  of  angn i t ode 


18 


3*-  '  <  *  *  -  •«. 


as  the  imposed  field. 

Linear  Regime 

Regardless  of  the  initial  write  beam  intensity,  the 
PRIZ  device  always  starts  in  the  linear  regime.  From 
equation  2.15,  the  excitation  rate,  g(0),  will  determine  the 
length  of  time  that  the  device  will  spend  in  the  linear 
regime.  The  linear  regime  is  characterized  by  an  electric 
field  which  can  be  approximated  everywhere  by  Bq  and  a 
stationary  interface  point.  Approximating  E  by  EQ  means 
that  E  is  no  longer  dependent  on  the  conduction  band 
density,  n,  so  equation  2.9  is  a  linear  equation,  thus  this 
regime  is  referred  to  as  linear. 

Since  the  longitudinal  electric  field  can  be 
approximated  everywhere  by  EQ ,  no  bottleneck  arises.  If  we 
define  the  total  charge  number  density  as 

p=n+,-n-n_  (2.17) 

and  replace  E  by  EQ,  then  the  charge  number  density  can  be 
represented  analytically  by  (8:687) 


gt 


o  - 


1  1  -  1  ) 


exp 


-z 


lEolTt 


-  7  exp 


[  ) 


(2.18) 
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where  7  =  cu<|E0|rt.  Since  by  definition,  the  charge  density 
i9  zero  at  the  interface,  setting  equation  2.18  to  zero  and 
solving  for  z  gives  the  location,  zQ,  of  the  interface. 
(6:1688) 


z 


o 


t  In (7) 

(  7  -  1  ) 


(2.19) 


where  L  is  the  drift  length,  </|E0|rt.  The  interface  zQ  is 
stationary  throughout  the  linear  time  regime. 

We  can  now  determine  the  electric  field  arising  from 
this  charge  distribution.  In  a  one-dimensional 
approximation,  the  normalized  longitudinal  field  can  be 
calculated  in  terms  of  the  normalized  charge  density  as: 


96  p  (z) 

dz  nQ 

where  8  -  |E|/|E0|.  Integration  yields 


(2.20) 


*<0)t 

8  s  -  — —  exp 

nQ(l-7) 

where  Cj  is  the  constant  of  integration.  For  a  given 
potential  constant  across  the  PRIZ,  the  value  of  Cj  is 
established.  This  allows  us  to  express  the  total 


•z 

t 


g(0)t 
nQ( 1-7) 


exp 


(-az) 


+  C- 


(2.21) 
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longitudinal  field  as 


where  dQ  is  the  length  of  the  crystal.  This  represents  the 
complete  solution  for  the  electric  field  in  the  linear 
regime . 

Linear  Regime .  Electric  Field  Dependence  on  the  Absorption 
Coefficient 

The  form  of  the  electric  field,  equation  2.22,  can  be 
simplified  by  considering  two  cases:  when  the  absorption 
coefficient  is  small,  and  when  the  absorption  coefficient  is 
large.  In  the  case  of  a  small  absorption  coefficient  which 
corresponds  to  uniform  illumination  and  g(z)  -  g(0),  Bryskin 
has  approximated  the  charge  distribution  given  by  equation 
2.18  (7:31) 


(2.23) 
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where  e(x)  is  the  Heaviside  function. 


Using  this 


approximation  for  o  and  following  the  same  procedures  (Hr 
the  complete  solution  for  o  ,  the  electric  field  is  given  by 


I  E  | 


gt ( zQ  -  z  ) 
nO*Tt 


*  l«„l 


2dn„t 
o  o 


z  £  z„ 
o 


IB 


2donoL 


(2.24) 


z  £ 


z 
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where  zQ  is  as  defined  by  equation  2.19.  Examining  equation 
2.24,  we  can  see  that  the  field  at  the  illuminated  face  of 
the  crystal  grows  linearly  with  time,  and  beyond  the 
interface,  zQ,  the  field  is  independent  of  position.  This 
holds  true  only  for  the  case  where  o  is  small. 

Now  consider  the  case  where  the  device  is  still 
operating  in  the  linear  regime  but  the  absorption 
coefficient  is  large.  There  is  no  equivalent  simplifying 
approximation  to  make  for  the  space  charge  density  as  there 
was  when  the  absorption  coefficient  was  small.  We  instead 
will  examine  equation  2.21  and  see  how  we  can  simplify  that 
expression.  In  cases  where  a  is  large,  we  can  discard  the 
term  exp(-adQ).  Therefore  we  have  the  following 
approximation  for  the  field 
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|E|  =  |  Eq  |  ♦  |B0|B<^exp(-z/4)  -  exp(-az)} 


|Bo,a/i-[«p(.do/il  -  i]  *  ^ 


(2.25) 


with 


g(0)t 

B  =  -  (2.26) 

“o'1”' 


Thus  we  have  an  approximation  to  what  the  field  will  look 
like  in  the  case  of  strong  absorption. 

Nonl inear  Regime 

In  the  case  where  the  total  charge  created  is  large 
enough  to  produce  a  self  field  on  the  order  of  the  applied 
field,  equation  2.9  becomes  non-linear  in  the  charged 
particle  density.  The  electric  field  increases  in  the 
vicinity  of  the  illuminated  electrode  and  obtains  a  minimum 
value,  less  than  EQ ,  at  the  interface  location,  zQ .  This 
minimum  leads  to  a  bottleneck  in  the  flow  cf  electrons.  At 
the  bottleneck  the  electron  velocity  is  severely  restricted 
and  the  relative  probability  of  being  trapped  in  this  region 
is  in  turn  raised.  This  increased  trapping  produces  a 
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migration  of  the  interface  toward  the  illuminated  electrode. 

The  rate  of  accumulation  of  the  integrated  longitudinal 
charge  in  the  crystal,  Q,  is  the  prime  driver  behind  the 
change  in  the  electric  field  value.  Q  is  actually  a  surface 
charge  density  generated  by  the  charge  distribution 
function,  p.  An  analytic  description  due  to  Bryskin  starts 
by  defining  how  the  total  charge  surface  changes  with  time 
(7:30) 


dQ 

—  =  j(0,t)  -  j (d  , t )  (2.27) 

dt  0 

where  j(0,t)  is  the  injection  current  and  j(d  ,t)  is  the 
current  leaving  through  the  dark  face.  From  equation  2.4 
the  current  j  can  be  written  as 

j  =  e«nE  (2.28) 

Bryskin  then  assumes  that  we  are  again  operating  where  the 
absorption  coefficient,  o,  is  small.  Using  the 
approximation  in  equation  2.23  to  find  the  surface  charge 
density,  Q,  yields 


Q  =  *tzQ 


(2.29) 


If  we  assume  the  injection  current  vanishes,  all  we  need 
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worry  about  is  the  current  traversing  the  dark  face.  We  can 
write  the  field  in  the  region  beyond  the  interface  using 
equation  2.24.  Substituting  this  field  into  equation  2.28 
and  then  substituting  equation  2.28  and  2.29  into  2.27 
yields 


d  r  1 

r  o  'i 

gtz  * 

=  *n(d0,t)|E0| 

1 - 2 - 

dt1-  °J 

2d„n. L 

L  o 

We  now  have  the  increase  in  the  surface  charge  density  in 
the  crystal  defined  by  the  current  leaving  the  crystal  in 
terms  of  the  electric  field  approximation  in  the  case  of  a 
small . 

Bryskin  next  finds  an  expression  for  the  electron 
conduction  band  density  at  the  dark  face.  Approximating 
equation  2.17  by 


p  =  gt  -  n_  -  n  (2.31) 

and  using  the  fact  that  in  our  approximation  p  -  0  beyond 
zQ,  allows  us  to  rewrite  equation  2.8  as 

4n  n 

—  +  —  =g  z  ^  zn  (2.32) 

at  t  t 

using  the  initial  conditions  n(zrO)  -  0,  the  solution  to 
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equation  2.32  is 


n  s 


(2.33) 


In  the  nonlinear  regime  t  >>  so  we  will  neglect  the 

exponential  term  in  equation  2.33.  Now  we  can  write 
equation  2.30  as 


d 

dt' 


f  1 

r  2  4 

*tzo 

tzJ  = 

1 - 2 - 

1  °J 

2d0nQ^ 

(2.34) 


At  large  times  the  solution  can  be  written  (7:32) 
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2V€  ' 


1/2 


[eg(0)t 


(2.35) 


Solving  equation  2.24  for  the  illuminated  face  of  the 
crystal  and  approximating  for  times  t  >>  tQ  yields  an 
expression  for  the  field  at  the  illuminated  face  of  the 
crystal  (7:32) 
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where  V  is  the  applied  potential,  and  €’  is  the  dielectric 
constant  for  the  crystal.  For  cases  when  a  is  large,  this 
analysis  is  inappropriate  and  numerical  methods  must  be 
used. 

Aa  cratov  Solution 

Now  we  have  identified  a  need  for  an  iterative 
numerical  solution  for  certain  cases  of  operation  of  the 
PRIZ,  notably  the  transition  region  between  the  linear 
regime  and  the  nonlinear  regime  and  the  nonlinear  regime  for 
cases  where  the  absorption  coefficient  is  large.  Astratov 
has  presented  a  numerical  solution  of  the  conduction  band 
density,  longitudinal  electric  field  and  space  charge 
density  assuming  a  quasi-steady  state  (3:1586): 

dn 

—  -  0  (2.37) 

dt 

Unlike  the  previous  analysis  of  Bryskin,  he  assumes  that  the 
source  term  g(z)  results  solely  from  thermal  excitation  of 
electrons  from  previously  uniformly  populated  shallow  traps. 
The  excitation  term,  g(z,t),  decays  with  time  as  the 
population  relaxes  thermally  with  relaxation  time, 
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Thus  (3:1585) 


g(t) 


g( 0)exp 


(2.38) 


With  the  assumption  made  in  equation  2.37  it  becomes 
possible  to  solve  equation  2.9  analytically  with  the 
boundary  conditions  specified  by  equations  2.12  and  2.13. 

The  solution  becomes  discrete  as  the  partial  derivative  with 
respect  to  z  is  approximated  by  Az.  Equation  2.9  is  set 
equal  to  zero  and  solved  by  employing  the  integrating  factor 
approach  (4:350).  Astratov  obtains  an  expression  for  n(z) 
at  fixed  points  in  the  crystal,  the  points  designated  by  the 
integer  i 


E(z,  , ) 

n(zi)  s  ■  ■n(zi,1) 

1  E(z.)  11 


exp 


Az 


s|E(zi) | r  t 


g(t>rt 


1  -  exp 


A  Z 


*(|E(zi)  |rt 


(2.39) 


where  E(z^)  is  the  field  at  the  specified  point  and  Az  is 
the  distance  in  between  specified  points.  The  total  charge 
number  density  is  found  by  adding  the  existing  density  at 
time  t  to  the  conduction  band  density  at  t+At  and  the  fixed 
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charge  at  t+At 


Ap  f  (  z i  ) 


n(Z;  ) 

g(t) - L- 


At 


(2.40) 


Therefore,  the  total  charge  number  den9ity  can  be  found  by 
the  following  expression 

p(zi,t+At)  =  p  (  z  ^  ,  t )  +  Apf(zi)  -  n(zi)  (2.41) 

The  electric  field  at  any  point  i  is  calculated  by  taking 
the  charge  density  at  all  points  i,  assuming  the  charge  to 
be  distributed  in  an  infinite  plane,  and  summing  each  planes 
contribution  to  the  electric  field  at  the  point  i.  Once  the 
new  electric  field  is  found,  the  g(t)  term  is  advanced  by  At 
and  the  solution  can  proceed  iteratively  through  time. 

The  Astratov  solution  is  flexible  enough  to  accept 
other  types  of  time  dependent  excitation  merely  by  rewriting 
g(t)  appropriately.  It  is  possible  to  extend  Astratovs' 
basic  solution  to  a  more  general  form  which  treats 
absorption.  In  order  to  treat  absorption,  it  is  necessary 
to  start  again  with  the  basic  set  of  equations  and  solve 
using  the  new  form  of  g(z).  Writing  the  excitation  rate  to 
include  absorption,  g(z)  =  g(0)  exp(-az),  the  solution  to 
n(z^)  is  found  to  be 
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,  ,  B(z*_.,) 

n(z*}  =  —  ■  *  *  n(2j  .)  exp 
1  E(Zi)  1-1 


Az 


it  |  B(  zi  )  |  rt 


g(0,t)rt  i  ( 

- - — *»— Kexp  -zta  - 

4 )rfca  +1]  1  1  J 


«E<2i>'t 


exp 


Az 


I  T  t 


“zi-la 


(2.42) 


The  solution  to  the  modified  Astratov  proceeds  as  outlined 
previously. 

The  Astratov  solution  is  not  valid  for  times  in  which 
the  conduction  band  density  is  rapidly  changing.  In 
addition,  the  method  provides  no  guidance  for  numerical 
stability,  the  solution  can  be  overdriven  if  the  time  step 
is  too  large,  and  no  provision  exists  in  either  the  Astratov 
solution  or  the  modified  Astratov  solution  to  correct  a  time 
step  that  is  too  large.  The  Astratov  solution  is  also 
limited  to  the  simplified  system,  equations  2.6-8,  if  an 
additional  term  is  to  be  examined,  the  system  must  be 
resolved,  if  indeed,  that  is  even  possible. 
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Chapter  III 


The  previous  chapter  identified  two  operating  regimes 
of  the  PHIZ.  This  chapter  presents  a  numerical  solution  of 
the  system  of  equations  2.7-9  which  is  applicable  in  either 
operating  regime.  The  method  used  to  solve  the  system  of 
equations  will  be  presented,  followed  by  comparisons  of  the 
numerical  solutions  with  analytic  results.  Normal  PR  I  Z 
operating  conditions,  however,  do  not  satisfy  many  of  the 
assumptions  on  which  these  analytic  solutions  were  based. 
Therefore,  the  numerical  model  will  be  used  to  examine 
solutions  under  these  conditions. 

Numerical  Mode  I  Deve 1 opment 

From  the  system  2.7-9,  we  can  develop  a  system  of  two 
equations  which  can  be  used  to  find  the  total  charge  in  the 
sample.  The  first  equation  of  interest  is  equation  2.9. 

4n  n  9  ( 

—  =  g(z)  -  —  ♦  a—  En  (3.1) 

9 1  9  Z  '  ' 

By  defining  a  new  term,  p^  :  nf  -  n_,  we  can  find  the  time 
rate  of  change  of  p^  by  subtracting  equation  2.7  from 
equation  2.8.  This  yields 


31 


g  (  Z  t 


n 


i  .  ^ 


r 


t 


Equation  3.1  contains  differentials  of  two  different 
variables,  z  and  t.  The  derivative  with  respect  to  z  *n  be 
approxiaated  as  a  difference  to  give  the  following 


in 


J  t. 


n 


,  t 


) 


n  [2  i  -  1  ’  t  ]  H  f  2  i  -  1 


l  3  •  3  1 


where  ~z  is  the  distance  between  the  discrete  bins 
designated  by  the  subscript  i.  Referring  back  to  Figure 
1.1,  we  have  divided  the  crystal  longitudinally  into  m 
slices,  where  the  index  nuabered  0  is  at  the  write  beam  side 
and  the  index  nuabered  a  is  at  the  opposite  side.  The  first 
tera  to  the  right  of  the  equals  sign  in  equation  3.3  is  the 
local  rate  of  excitation  into  the  conduction  band.  The 
second  tera  is  the  rate  of  electron  decay  out  of  the 
conduction  band  due  to  trapping.  The  last  tera  in  the 
brackets  approxiaates  the  divergence  of  the  conduction  band 
electron  flux.  We  thus  have  a  systea  of  a  coupled, 
first-order,  nonlinear  equations  which  describe  the  density 
of  conduction  band  electrons  in  the  crystal.  We  rewrite 
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equat ion  3.2  as 


Here,  g(z^)  represents  the  rate  of  production  of  fixed 
positive  charge  sites  at  the  position  z^.  The  second  term 
represents  the  formation  rate  of  fixed  negative  charge 
sites.  The  total  net  charge  at  any  given  time  can  be  found 
by  subtracting  the  conduction  band  density  from  the  net 
fixed  charge. 


(3.5) 


Once  the  charge  density  at  every  point  is  known,  the 
electric  field  at  every  point  can  be  calculated.  This  is 
most  easily  done  in  a  one  dimensional  approach  by  assuming 
the  charge  density  at  each  point  is  distributed  in  an 
infinite  plane.  The  value  of  the  electric  field  at  any  point 
z^f  is  found  by 


E 


(Z 


i 


i 


m 


k  =  1 


(3.6) 
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where  a 


jj  =  fl  (Z|j)Ar,  the  area  charge  density  at  any  discrete 
point  k,  Ar  is  the  area  of  the  PRIZ  crystal,  £  is  the 
dimensionless  dielectric  constant  of  the  crystal  and  £  is 
the  dielectric  constant  of  free  space.  The  charge  density 
and  electric  field  are  calculated  self-conaistently  at  each 
time,  t. 

The  IMSL  Fortran  subroutine,  DGEAR  (11:8),  was  chosen 
to  solve  the  equations  3.3  and  3.4  as  a  system  of  first 
order  differential  equations.  DGEAR  was  appropriate  because 
the  subroutine  uses  backward  differentiation  for  stiff 
equations.  Also,  DGEAR  is  designed  to  handle  an  initial 
condition  problem,  and  the  initial  conditions  of  the  PRIZ 
are  very  easy  to  write, 


n[zt,o]  =  0 


(3.7) 


P 


=  0 


(3.8) 


The  spatial  boundary  conditions  of  the 
were  specified  by  modifying  equation  3.3  at 
i  s  l  and  ism.  The  boundary  condition  at 
face  was  specified  by  eliminating  the  flux 


conducting  PRIZ 
the  points  where 
the  illuminated 
into  the  first 


bin,  equation  2.12.  The  boundary  condition  at  the  dark  face 


was  specified  by  allowing  the  free  escape  of  conduction  band 
electrons  into  a  dummy  discrete  point,  m+1. 

Conduct ing  PRIZ  Device 

The  device  modeled  in  the  remainder  of  the  chapter  is  a 
PRIZ  constructed  of  a  (111)  cut  BSO  crystal.  The  dielectric 
constant,  t,  is  56.0.  The  mobility  of  the  electrons,  u,  is 
3.0  x  10”®  m^/(sec-V).  The  lifetime  of  the  electrons  in  the 
conduction  band,  rfc,  is  3.0  X  10  sec.  The  crystal 
parameters  listed  above  are  taken  from  ref.  8.  The  1cm  * 
lcm  x  1mm  crystal  has  a  potential  of  2000V  applied  across 
the  1mm  dimension.  It  is  assumed  that  this  potential 

difference  is  maintained  throughout  the  period  of  operation. 

_  2 

This  yields  a  drift  distance,  L  =  1.8  x  10  cm.  Three 
write  beam  wavelengths  are  used  in  order  to  examine  the 
sensitivity  to  wavelength  variations  in  the  absorption 
coefficient.  Weak  absorption  corresponds  to,  x  =  589  nm  and 
a j  =  0.6  cm”* ;  strong  absorption,  X  =  400  nm  and  =  90.0 
cm-1.  Intermediate  comparisons  at,  X  =  434  nm  and  a „  =  45.0 
cs'*  are  also  included  (17:13).  The  data  presented  assumes 
that  the  excitation  rate  at  the  illuminated  face  of  the 
crystal,  g(0),  is  constant  for  ail  cases.  This  results  in 
the  excitation  of  a  larger  number  of  charge  carriers  as  the 
absorption  coefficient  becomes  smaller.  Conversions  back  to 
incident  power  can  be  made  using  equation  2.14. 
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L  i  near  Re  g  1  me  .  C2*D®rL52!ll-  with  Coaplet^e  Solutions 

For  the  Linear  regiae,  where  t  <<  tQ,  there  exist 
analytic  foraa  which  can  be  used  to  predict  both  the  spa*  »? 
charge  density,  equation  2.18,  and  the  longitudinal  electric 
field,  equation  2.22.  For  caaea  in  which  the  absorptiwn 
coefficient,  a,  ia  aaall,  theae  equations  simplified  to 
charge  density,  equation  2.23,  and  electric  field,  equation 
2.24,  respectively.  When  the  absorption  coefficient  is 
large,  the  electric  field  was  shown  to  be  approximated  by 
equation  2.25.  Comparisons  of  the  numerical  results  with 
the  analytic  forms,  and  both  sets  of  approximations  will  be 
presented  in  this  section.  In  order  to  insure  that  the 
linear  regime  approximation  is  valid,  t  -  .01to>  for  all 

compar i sons . 

The  first  comparison  to  consider  is  that  of  the 
numerical  solution  to  the  exact  analytic  forms  which  exist. 
Figures  3.1-3  compare  the  predictions  of  equation  2.1H  to 
the  results  yielded  by  the  numerical  solution  for  small, 
medium  and  large  absorption  coefficients.  The  curves  .an  be 
seen  to  exhibit  excellent  agreement,  the  point  of  greatest 
disparity  is  in  the  prediction  of  the  point  zQ.  For  Uj  the 
difference  as  a  fraction  of  the  total  crystal  length  is 
5.36%.  For  a^,  the  difference  is  1.57%.  For  a ^ ,  the 
difference  is  1.57%. 

The  surface  charge  number  density  of  the  PRI2,  ij,  is 
the  determining  factor  in  the  strength  of  the  self  electric 
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Figure  3. 


1.  Linear  Regiae,  Longitudinal  Charge  Number 
Density  for  small  Absorption  Coefficient, 
o  s  0 . 6  ca“  ,  t  s  . 0 1  tQ 
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Charge  Number  Density 


Figure  3.2.  Linear  Regime,  Longitudinal  Charge  Number 
Density  for  medium  Absorption  Coefficient, 
a  =  45  cm"1,  t  =  .01tQ 
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Figure  3. 


3.  Linear  Regime,  Longitudinal  Charge  Number 
Density  for  large  Absorption  Coefficient, 
a  =  90  cm"  ,  t  =  . 01tQ 
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field.  To  analytically  calculate  the  surface  charge  number 
density,  Q,  we  integrate  equation  2.18,  which  yields 

g(0)t L 

Q  =  - 

1  -  7 

From  the  numerical  solution,  the  surface  charge  number 
density  can  be  numerically  integrated  using  the  trapezoid 
rule.  For  a^,  the  difference  between  the  analytic  form  and 
the  numerical  Q  is  0.24%.  The  analytic  value  for  Q  is 
slightly  larger  than  the  numerical  value  for  Q.  For  a3 ,  the 
difference  is  2.91%,  with  the  larger  value  being  predicted 
by  the  numerical  solution.  For  a3,  the  difference  is  8.3%. 
The  numerical  solution  again  predicts  the  larger  value. 

These  results  agree  quite  well  and  are  expected  from  the 
correlation  of  the  charge  density  curves. 

Now  we  will  compare  the  predictions  of  the  complete 
solution  to  the  longitudinal  electric  field,  equation  2.22, 
with  the  numerical  solution.  Equation  2.22  is  derived  using 
the  analytic  form  of  the  charge  density,  equation  2.18.  The 
total  field  E  is  not  compared,  rather,  the  self  field 
generated  by  the  internal  charge  distribution  only  will  be 
plotted.  The  self  field  at  the  illuminated  face  of  the 
crystal  can  be  calculated  directly  from  the  total  charge 
number  density,  Q.  The  greatest  disparity  in  the  solutions 


exp[-adQ]  -  exp[-dD/^j 


3.9) 
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occurs  in  calculating  the  self  field  at  the  illuminated  face 
of  the  crystal.  For  a^,  the  numerical  solution  is  the 
larger  of  the  two  and  the  difference  between  them  is  18.9%. 
For  i  the  numerical  solution  is  again  larger  and  the 
difference  is  7.0%.  For  Og,  the  numerical  solution  is 
larger  and  the  difference  is  17.9%. 

The  numerical  solution  is  sensitive  to  the  grid  size 
chosen.  In  order  to  demonstrate  this,  the  number  of  bins 
was  increased  to  2m,  in  the  case  where  a  =  ctj.  This 
resulted  in  a  tenfold  increase  of  computer  time.  The 
agreement  between  the  values  for  Q  is  0.15%.  For  the 
electric  field  at  the  illuminated  face  of  the  crystal  the 
difference  is  0.96%.  Figure  3.4  displays  the  numerical 
solution  for  the  electric  field  only,  since  the  two 
solutions  are  essentially  the  same.  Figure  3.5  and  3.6  show 
the  comparison  between  the  numerical  forms  and  the  form  of 
equation  2.22.  Both  Figure  3.5  and  3.6  were  calculated  with 
m  discrete  points. 

Figure  3.4  shows  the  result  when  the  crystal  contains 
an  essentially  positive  charge  distribution  is  present 
throughout  the  crystal.  The  charge  distribution  for  this  is 
shown  in  Figure  3.1.  Figures  3.5  and  3.6  show  the  results 
for  the  crystal  when  there  is  a  positive  region  and  a 
negative  region,  these  charge  distributions  are  shown  in 
Figures  3.2  and  3.3,  respectively.  The  field  actually  forms 
a  slight  bottleneck  and  then  increases  again  due  to  the 
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Linear  Regime,  Longitudinal  Self  Electric 
Field  for  o  s  0.6  cm-1,  t  =  .01tQ 


Figure  3.5.  Linear  Regime,  Longitudinal  Self  Blectric 
Field  for  o  a  45  cm-1,  t  =  .01to 
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Figure  3. 


6.  Linear  Regime,  Longitudinal  Self  Electric 
Field  for  a  =  90  cm.  ,  t  =  .01to 


44 


negative  charge  region. 

The  analytic  forms  that  Bryskin  derives  for  the  linear 
regime  work  quite  well  when  the  time  under  consideration  13 
well  within  the  linear  regime.  The  complete  solution  for 
the  electric  field,  equation  2.22,  derived  in  this  study 
matches  excellently  with  the  numerical  solution.  The 
numerical  solution  produces  acceptable  results  in  a  short 
amount  of  computer  time.  If  computer  time  is  not 
restricted,  the  results  available  to  the  numerical  solution 
are  excellent. 

Linear  Regime .  Comparisons  with  Approximate  3o lut i ons 

The  analytic  solutions  developed  in  Chapter  II  were 
analyzed  for  their  behavior  in  limiting  cases  of  the 
absorption  coefficient.  These  limiting  cases  were  when  a  is 
small  and  when  a  is  large.  The  numerical  solution  developed 
in  this  chapter  will  be  compared  to  the  limiting  forms. 

The  numerical  solution  can  be  used  to  find  the  space 
charge  density,  o{ z),  in  the  limiting  case  of  a  small.  The 
approximation  for  p(z)  is  given  by  equation  2.23.  L'sing  a  j 
as  the  absorption  coefficient  that  is  small,  Figure  3.7 
shows  the  curves  generated  by  both  the  numerical  solution 
and  the  approximation.  The  form  of  these  curves  does  not 
compare  favorably,  so  the  best  way  to  check  for  the  validity 
of  the  approximation  is  to  compare  the  value  of  the  total 
charge  on  the  PRIZ,  Q.  Equation  2.23  can  be  seen  to  predict 
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a  larger  value  of  Q  than  that  of  the  numerical  solution. 

The  approximate  value  of  Q  ia  4.8  times  the  value  of  Q 
yielded  by  the  numerical  solution.  This  does  not  agree 
nearly  as  well  as  the  analytic  expression  given  in  equation 
2.18  does . 

Using  the  charge  density  approximation,  equation  2.23, 
the  electric  field  can  be  calculated.  This  expression  is 
given  in  equation  2.24.  A  comparison  of  the  two  curves  is 
shown  in  Figure  3.8.  The  value  of  the  electric  field  at  the 
illuminated  face  of  the  crystal  is  clearly  larger  for  the 
approximation  than  that  of  the  numerical  solution.  This  is 
expected  due  to  the  fact  that  the  approximation  severely 
over  predicted  the  total  charge  present  in  the  crystal.  The 
value  of  the  approximate  solution  is  1.72  times  the 
numerical  solution. 

Clearly,  this  is  not  a  good  approximation  for  even 
times  well  within  the  linear  regime.  Recall  that  we  are 
operating  at  .01tQ,  where  tQ  is  the  temporal  boundary 
between  the  linear  and  nonlinear  regimes.  Recalling  the 
development  of  the  nonlinear  equations  at  the  end  of  chapter 
II,  the  motivation  for  the  approximation  becomes  clearer. 
Bryskin  desired  a  form  that  would  lend  itself  to  analytic 
solution  when  determining  the  forms  of  th>  electric  field  at 
the  illuminated  face  of  the  crystal  and  the  location  of  the 
interface,  zQ. 

For  the  case  where  the  absorption  coefficient  is  large, 
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Figure  3.8.  Linear  Regine,  Longitudinal  Self  Electric 
Field  for  a  =  0.6  cm*  ,  t  =  .01tQ 
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no  approximation  for  the  apace  charge  density  is  given.  The 
comparison  between  the  analytic  fora  of  the  space  charge 
density  and  the  numerical  solution  of  space  charge  density 
is  shown  in  Figure  3.3. 

The  approximate  form  of  the  electric  field  in  the  case 
where  a  is  large  is  given  in  equation  2.25.  Since  equation 
2.25  is  very  similar  to  equation  2.22,  the  difference 
between  the  value  of  the  electric  field  at  the  illuminated 
face  of  the  crystal  for  the  approximate  solution  and  the 
numerical  solution  should  be  very  close  to  the  difference 
found  between  the  numerical  and  analytic  cases.  For  large 
a,  using  as  large  a,  the  difference  is  18. OX.  Figure  3.9 
shows  the  numerical  and  approximate  solution  for  a  large 
absorption  coefficient. 

Nonlinear  Regime 

For  times  greater  than  tQ,  we  can  no  longer  approximate 
the  electric  field  by  EQ.  The  nonlinear  section  in  Chapter 
II  developed  solutions  for  the  electric  field  at  the 
illuminated  face  of  the  crystal  and  the  location  of  the 
interface,  zQ,  under  the  assumption  of  a  small  absorption 
coefficient.  Therefore,  we  would  not  expect  these  results 
to  hold  for  large  absorption  coefficients.  It  is  also 
useful  to  remember  that  the  approximation  used  for  the 
charge  number  density  2.23  over  predicted  the  total  Q  of  the 
crystal  by  a  factor  of  4.8.  So  even  in  the  case  where  the 
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Fi*ure  3.9.  Linear  Regime,  Longitudinal  Self  Electric 
Field  for  a  s  90  c«  ,  t  s  .01tQ 
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absorption  coefficient  is  small,  the  prediction  of  the 
electric  field  at  the  illuminated  face  and  the  location  of 
the  charge  interface  could  be  seriously  in  error.  The 
predictions  are  more  useful  for  the  predictions  of  the  time 
dependences  of  the  parameters. 

The  electric  field  at  the  illuminated  face  of  the 

crystal  for  both  the  approximate  and  the  numerical  solution, 

is  shown  in  Figure  3.10  for  times  t  >  1.5tQ,  and  where  a  is 

small,  =  0.6.  The  results  of  equation  2.36  are 

approximately  a  factor  of  1.5  greater  than  the  numerical 

results.  This  is  better  than  we  anticipated  in  view  of  the 

comparison  of  predicted  surface  charge  number,  Q.  The 

0  7  4 

temporal  correlation  exhibits  a  dependence  of  t  *  ,  in 

contrast  with  the  analytic  result  of  t0,5. 

Figure  3.11  shows  the  electric  field  at  the  illuminated 
face  of  the  crystal  for  both  the  approximate  and  the 
numerical  solution  for  the  same  time  as  Figure  3.10  but  with 
otg  =  90.0.  Notice  that  equation  2.36  does  not  depend  on  the 
absorption  coefficient.  Thus  the  results  of  equation  2.36 
are  approximately  a  factor  of  7.0  larger  than  the  numerical 
solution.  The  time  dependence  of  the  numerical  solution  is 
now  The  fit  for  both  the  magnitude  of  the  field  and 

the  time  dependence  of  the  field  exhibit  poor  agreement. 

This  is  to  be  expected  because  1 )  in  all  cases  the 
approximation  from  which  the  nonlinear  results  derived, 
equation  2.23,  overestimates  the  charge  present  in  the 
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Figure  3. 


0.  Nonlinear  Regime,  Magnitude  of  |E(0,t)| 
with  respect  to  time  for  a  s  0.6  cm'1 
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Figure  3.11.  Nonlinear  Regime*  Magnitude  of  |E(0.t)| 
with  respect  to  time  for  a  =  90  cm 
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crystal;  and  2)  the  absorption  coefficient  was  assumed  to  be 
small  throughout  the  derivation  of  the  nonlinear  results. 

The  PRIZ  write  beamf  however,  is  usually  of  such  a 
wavelength  so  as  to  cause  the  absorption  coefficient  to  be 
large.  In  addition,  the  beam  is  usually  intense  enough  to 
cause  the  device  to  quickly  enter  the  nonlinear  regime. 

Thus,  for  most  real  applications  of  the  device,  only  the 
numerical  solution  will  suffice. 

The  approximate,  equation  2.35,  and  numerical  solutions 
for  the  location  of  zQ  at  various  times  when  the  absorption 
coefficient  is  small,  ,  are  presented  in  Figure  3.12.  The 
numerical  results  exhibit  a  dependence  of  t“*555,  which 
compares  favorably  with  the  analytically  predicted  t-0’® 
dependence.  The  analytic  expression  predicts  the  value  of 
zQ  to.be  approximately  0.8  of  the  numerical  solution. 

A  similar  comparison,  now  with  a  large  absorption 
coefficient,  ag,  is  presented  in  Figure  3.13.  The  numerical 
data  yields  a  dependence  of  for  t  >  1.5t0. 

Bryskins’  value  of  zQ  is  roughly  2.0  times  as  far  away  from 
the  illuminated  face  as  the  numerical  prediction.  The 
disparity  in  the  location  of  zQ  is  expected,  since  Bryskins’ 
derivation  assumed  that  the  absorption  coefficient  was 
small.  The  agreement  in  the  time  dependence  of  zQ  is 
somewhat  surprising. 

We  have  seen  that  when  t  <<  tQ ,  the  analytic 
expressions  for  the  charge  number  density,  equation  2.18, 
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Figure  3.12.  Nonlinear  Regime,  Location  of  Charge 

Interface,  zQ,  with  respect  to  time  for 
o  =  0.6  cm"1 
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Distance  of  z  from  Crystal  Face  (ca  xlO 


3.34 


Figure  3.13.  Nonlinear  Regime,  Location  of  Charge 

Interface,  z  ,  with  respect  to  time  for 
o  s  90  cm’1 
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the  surface  charge  number  density,  equation  3.9,  and  the 
electric  field,  equation  2.22,  exhibit  excellent  agreement 
with  the  numerical  results.  This  agreement  holds  regardless 
of  the  absorption  coefficient  of  the  material  for  the 
wavelength  under  consideration.  When  t  >  tQ,  the  analytic 
expression  for  the  charge  interface,  equation  2.35,  exhibits 
acceptable  agreement  for  all  absorption  coefficients.  The 
expression  for  the  electric  field  at  the  face  of  the 
crystal,  equation  2.36,  is  only  useful  when  the  absorption 
coefficient  is  small. 
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Chapter  IV 


In  the  previous  chapter  the  numerical  model  of  the  PRIZ 
was  developed  and  validated  using  analytical  results.  It  is 
■ore  difficult  to  evaluate  the  model  when  considering 
experimental  results.  Soviet  data  on  the  longitudinal 
electric  field  exists.  However,  these  data  omit  critical 
parameters  that  the  numerical  model  requires  for  input. 

’sin*  best  estimates  for  these  parameters,  a  comparison  is 
undertaken ,  with  any  assumptions  made  clearly  detailed. 

Two  interesting  experimental  observations  relating  to 
PRIZ  operating  characteristics  are  reported  by  Shields 
ref.  17)  and  by  Anderson  (ref.  1  and  ref.  2).  Shields 
reports  on  the  "edge  enhancement"  of  a  simple  spot.  The 
'■barge  distribution  provided  by  the  numerical  model  are  used 
*■>  calculate  the  transverse  fields  in  this  simple  case,  and 
*  he  result  is  related  to  Shields’  observations.  Anderson  has 
reported  that  the  BSO  crystal  suffers  damage  when  operating 
in  the  dynamic  imaging  mode.  Observations  from  the  model 
will  be  offered  to  help  explain  this  phenomenon. 

. ompar  t son  o  f  Model  with  Experiment 

Bryskin  offers  experimental  results  in  two  papers, 
ref.  6  and  ref.  8).  These  results  take  on  the  form  of  the 
spatial  dependence  of  the  longitudinal  field  for  a  given 
set  of  parameters.  Unfortunately,  Bryskin  does  not 
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completely  specify  his  operating  conditions. 

The  results  reported  in  ref.  6  are  compared  with  the 
numerical  model  in  Figure  4.1.  The  total  exposure  was  not 
reported,  nor  was  the  time  at  which  the  measurements  were 
taken.  However,  it  can  be  assumed  that  as  long  as  the  same 
total  energy  is  deposited  on  the  PR1Z,  the  write  time  should 
be  immaterial  to  the  shape  of  the  longitudinal  electric 

C 

field.  Curve  1  has  an  applied  field  of  1.3  x  10  V/m  and 

g 

curve  2  has  an  applied  field  of  3.2  x  10  V/m.  The  power  of 

2 

the  write  beam  at  441  nm  was  assumed  to  be  1.11  sW/cm  and 

the  time  of  best  fit  was  observed  to  be  4.0  seconds,  for  a 

2 

deposition  of  4.44  uJ/cu  .  The  greatest  disparity  can  be 
seen  in  the  electric  field  behavior  beyond  the  bottleneck. 

For  curve  1,  the  numerical  result  predicts  a  slight  rise  in 
the  magnitude  of  the  electric  field,  while  the  experimental 
data  do  not  reflect  this.  For  curve  2,  the  rise  predicted 
by  the  numerical  result  is  matched  by  the  experimental  data, 
although  the  magnitude  of  the  rise  for  the  experimental  data 
is  less  than  the  rise  in  the  numerical  data. 

Figure  4.2  shows  the  results  reported  in  ref.  8 
compared  with  the  numerical  model.  The  write  intensity  was 
again  not  specified  in  this  model,  but  the  time  of  the 
measurement  was,  therefore,  once  a  reasonable  match  was 
achieved  for  a  given  excitation  rate,  g<0),  the  power  of  the 
write  beam  at  589  nm  could  be  calculated  to  be  374  *W/cn^, 
from  equation  2.14.  It  should  be  noted  that  unlike  in  Figure 
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i.  Longitudinal  Electric  Field,  Numerical  vs. 
Experimental  (6:1687) 
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4.1,  where  the  quantum  efficiency  was  1.0  and  the  absorption 
coefficient  41.0  cm” * ,  the  quantum  efficiency  in  Figure  4.2 
was  0.03  and  the  absorption  coefficient  0.6  cm“^.  So  while 
almost  all  of  the  energy  was  absorbed  in  the  experiment 
represented  by  Figure  4.1,  very  little  of  the  energy  was 
actually  aosorbed  by  the  crystal  in  the  experiment  in  Figure 

4.2.  In  Figure  4.2  the  agreement  between  experiment  and  the 
numerical  model  is  much  closer  throughout  the  crystal  than 
it  was  for  Figure  4.1.  The  possible  exception  is  the  value 
of  the  field  at  the  illuminated  face  of  the  crystal. 

For  these  two  cases,  the  numerical  model  was  easily 
able  to  match  the  experimental  results  to  within  a  factor  of 
two.  Thus  the  model  should  prove  valuable  as  a  tool  for 
experimenters  in  developing  further  experiments  dealing  with 
the  PRIZ  device. 

The  Case  of  a  Simple  Spot 

In  the  first  chapter,  an  expression  for  the  phase 
difference  of  the  PRIZ  crystal  was  developed,  equation  1.1. 
This  expression  was  dependent  on  the  integrated  transverse 
electric  field.  Up  to  now,  we  have  been  exploring  solutions 
for  the  longitudinal  space  charge  density  and  the 
longitudinal  electric  field.  It  was  necessary  to  consider 
these  together  since  they  represent  a  coupled  system.  Now 
that  we  know  the  space  charge  density,  we  may  now  find  the 
integrated  transverse  field  for  a  simple  case. 


62 


Method  of  Solution  f or  a  Simple  Spot 

The  case  to  be  treated  is  a  circular  write  beam  of 
small  radius  illuminating  a  spot  on  the  crystal  and 
generating  a  cylindrical  charge  distribution  within  the 
device.  From  this  cylinder  of  charge ,  equipotential 
surfaces  in  the  radial  direction  are  calculated.  The 
equipotential  surfaces  are  approximated  by  treating  the 
cylinder  of  charge  as  a  set  of  discrete  point  charges  and 
solving  for  the  potential  at  discrete  points  in  the  radial 
plane.  Thus  in  this  solution,  at  any  given  point  in  the 
potential  matrix,  the  potential  contributions  from  the  m 
discrete  points  in  the  longitudinal  cylinder  must  be  summed. 

Once  the  numerical  representation  of  the  potential  due 
to  the  charged  cylinder  is  found,  the  transverse  electric 
field  is  found  by  numerically  approximating  the  differential 
between  points  in  the  potential  matrix.  This  gives  us  a 
field 


E 


t 


dv  Aw  W  j  -  v 
97  '  Ay  "  -  7j _! 


(4.1) 


where  Et  is  the  transverse  electric  field  at  point  j,  and  w 
is  the  potential  at  the  given  point.  The  effect  of  the 
transverse  field  on  the  electrons  in  the  conduction  band  is 
ignored.  Further,  in  this  case,  the  write  beam  is 
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considered  to  have  a  very  small  radius ,  on  the  order  of  1 
um,  and  all  calculations  are  made  outside  the  radius  of  the 
beam.  Hence,  no  effort  was  made  to  model  the  transverse 
electric  fields  in  the  actual  write  beam  as  it  passes 
through  the  crystal.  Since  for  equation  1.1  we  need  to  know 
the  integrated  value  of  transverse  field,  the  transverse 
fields  at  every  given  radial  distance  y  from  the  write  beam 
are  integrated  along  the  longitudinal  axis. 

Results  in  the  Case  of  the  Simple  Spot 

Figure  4.3  is  a  representative  example  of  the 
equipotentials  due  to  the  write  beam  inside  the  crystal. 

The  closely  spaced  lines  show  the  concentration  of  positive 
charge  near  the  illuminated  face  of  the  crystal.  From  this 
equipotential  we  can  find  the  transverse  fields  at  the 
discrete  points  for  which  the  potentials  are  known.  Figure 
4.4  is  a  representative  example  of  the  transverse  fields  in 
the  crystal.  The  horizontal  axis  of  the  figure  labeled 
longitudinal  axis  is  the  direction  parallel  to  the  write 
beam.  The  vertical  axis  is  the  magnitude  and  direction  of 
the  transverse  field  at  that  radial  distance.  The  axis 
facing  more  or  less  out  of  the  page  is  the  radial  distance 
from  the  beam.  Figure  4.4  shows  that  the  largest  magnitude 
of  the  field  occurs  closest  to  the  beam  and  is  generated  by 
the  large  concentration  of  positive  charge  near  the 
illuminated  face  of  the  crystal. 
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Figure  4.3.  Equipotentiela  for  •  See 11  Cirouler 

Write  Beea 

•  ft 
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Radial  Longitudinal 

Axis  Axis 


Figure  4.4.  Transverse  Electric  Fields  for  a  Snail 
Circular  Write  Bean 
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In  order  to  find  the  integrated  value  of  the  transverse 
fields  in  the  longitudinal  direction,  each  of  the  curves 
representing  the  magnitude  and  direction  of  the  transverse 
field  at  a  fixed  radial  distance  along  the  longitudinal  axis 
of  beam  is  integrated.  Figure  4.5  displays  the  results  of 
this  integration  by  examining  the  integral  value  with 
respect  to  radial  distance  from  the  beam.  The  dotted  line 
present  in  the  graph  represents  a  strict  y-*  dependence.  To 
a  reasonable  approximation  then,  we  see  that  the  magnitude 
of  the  integrated  value  of  the  beam  has  an  approximately 
inverse  relationship  with  the  radial  distance  from  the  beam. 

Shields  reports  what  he  calls  the  "Edge-Enhanced  Output 
Image  of  Input  Spot"  (17:22).  The  results  shown  in  Figure 
4.5  would  indicate  that  we  would  find  the  largest  value  of 
the  integrated  transverse  field,  Bt’,  at  or  near  the  edge  of 
the  spot.  From  equation  1.1  the  value  of  E^.’  can  be  seen  to 
be  directly  proportional  to  the  phase  difference.  Since  rg^ 
is  always  a  very  small  number  the  phase  difference  will  most 
likely  never  exceed  w/2,  so  an  increase  in  the  value  of 
Et’would  bring  about  an  increase  in  the  intensity.  So  the 
edge-enhancement  is  something  we  would  expect. 

BSO  Crystal  Damage 

Anderson  reported  damage  occurring  to  the  BSO  crystal 
when  the  crystal  is  operated  in  the  dynamic  imaging  mode. 
O’Dwyer  (ref.  13)  reports  on  the  strength  of  electric  field 
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Figure  4 
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.5.  Integrated  Tranaverae  Electric  Field  for  a 
Saail  Circular  Write  Beam 
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necessary  to  cause  breakdown  in  various  substances.  Almost 

8 

all  of  his  results  fall  into  the  range  of  MV/cm  or  10  V/m. 
We  will  accept  this  figure  as  a  reasonable  order  of 
magnitude  estimate  for  the  field  required  to  cause  breakdown 
in  a  BSO  crystal.  Looking  back  at  Figure  4.1,  which 
corresponds  to  a  write  beam  at  wavelength  441  nm,  we  can 
compare  with  some  of  Andersons  results,  which  were  obtained 
with  a  442  nm  laser  source. 

2 

For  a  write  beam  of  1.11  *tW/cm  ,  after  4  seconds  the 
electric  field  at  the  illuminated  face  of  the  crystal  has 

C 

risen  to  4.7  x  10  V/m.  Increasing  the  write  beam  power 

density  by  a  factor  of  10  produces  an  electric  field  at  the 

illuminated  face  of  the  crystal  of  2.8  x  10^  V/m.  In  his 

initial  examination  of  damage  to  the  PRIZ  crystal,  Anderson 

uses  as  his  lowest  write  power  293  MW/cm  (1:4.5).  It  is 

2 

easy  to  see  that  a  write  power  of  293  sW/cm  could  produce  a 

a 

field  at  the  illuminated  face  in  the  range  of  10  V/m.  Thus 
it  is  highly  likely  that  the  observed  damage  resulted  from 
breakdown  in  the  vicinity  of  the  illuminated  surface. 
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Chapter  V 


Summary 

The  PRIZ  operates  by  creatine  a  phase  difference 
between  the  indices  of  refraction  perpendicular  to  the  axis 
of  propagation  of  the  write  and  read  beams.  This  phase 
difference  can  be  calculated  and  is  found  to  be  proportional 
to  the  strength  of  the  transverse  electric  field.  This 
field  arises  from  a  space  charge  generated  within  the 
crystal.  Thus,  establishing  the  distribution  of  the  space 
charge  is  a  principal  goal. 

This  study  presents  a  review  of  the  theory  governing 
the  excitation,  charge  transport  and  trapping  of  charge 
within  the  crystal.  The  theory  is  ail  baaed  on  the 
simplified  system  presented  in  equations  2.6-8.  A 
qualitative  discussion  of  the  processes  is  presented.  Two 
operating  regimes,  linear  and  nonlinear,  were  identified. 

The  theory  for  the  linear  regime  has  been  fairly  well 
established  by  Bryakin,  but  the  opportunity  was  taken  here 
to  do  more  rigorous  calculations  to  arrive  at  complete 
solutions  for  the  longitudinal  characteristics  of  the 
device.  The  numerical  method  used  in  solving  the  simplified 
system  of  coupled,  first-order,  differential  equations  is 
presented,  and  found  to  agree  quite  well  with  theory  in  the 
1 1  near  reg ime . 

Bryskin  has  also  developed  a  theory,  albeit  less 
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complete,  for  the  nonlinear  regime  of  the  device, 
theory  is  based  on  simplifying  assumptions  which  severely 
limit  the  quantitative  and  even  qualitative  utility  of  the 
results.  No  analytical  extensions  of  the  theory  in  this 
regime  was  attempted  in  this  study.  The  numerical  solution 
is  capable  of  operating  quit  well  in  either  regime. 

Experimental  results  reported  in  the  literature  are 
troubled  by  the  omission  of  several  important  parameters  in 
the  presentation  of  the  data.  The  experimental  data  which 
is  presented  is  compared  with  numerical  results,  the 
assumptions  made  in  the  model  in  order  to  complete  the 
parameter  set  are  clearly  laid  out.  Qualitative  results  are 
all  in  agreement.  Some  of  the  qualitative  observations 
reported  and  confirmed  by  the  numerical  solution  are: 

1.  The  longitudinal  electric  field  at  the  face  of 
the  crystal  increases  with  time  as  the  device 
is  written  upon. 

2.  The  longitudinal  electric  field  develops  a 
bottleneck  in  the  crystal  which  slows  the 
escape  of  charge  from  the  device. 

Shields  reports  the  edge  enhancement  of  a  simple  spot. 
The  numerical  solution  predicts  the  integrated  transverse 
field  at  the  edge  of  the  spot  to  be  a  maximum.  It  also 
predicts  the  essentially  1/r  fall  off  of  that  field. 

Anderson  (ref.  1  and  2)  has  studied  the  damage  to  the 
■  SO  crystal  when  operating  in  the  dynamic  imaging 
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Th**  postulated  mechanism  leading  to  the  damage  is 


<  •►••i  to  be  dielectric  breakdown.  Using  representative 

’"nu.’s  '»*en  from  O'Dwyer  (ref.  13),  the  numerical  model 
rhat  large  write  beam  powers  can  create  an  electric 
*  .« .  1  a ’  *  he  illuminated  face  of  the  crystal  which  is  large 

*>n  ugh  to  cause  breakdown. 

tr.  .mmendationa  for  Further  Research 

v. though  the  PRIZ  was  the  device  of  primary  interest 

*  study,  the  differential  equations  representing  the 
’ransport  in  the  BSO  crystal  are  general.  The  method 

*  .on  of  this  system  could  be  extended  to  other 
-  •  r**iymg  on  internal  charge  distributions  to  generate 

n<i tudinal  or  transverse  fields. 

i  9:  1  347)  has  proposed  using  the  PRIZ  for 
»»  mage  selection  and  Nilius  (ref.  12)  has  examined 

>-t  ’-m  memory  of  the  device.  Both  of  these  uses  depend 

•  m -  required  for  optical  erasure.  Therefore,  it  is 

*  •  ' mr  the  mechanism  causing  optical  erasure  be 

■<iiius  reported  the  decay  time  of  the  image, 

•  •  •  related  to  the  incident  power  on  the  crystal  by 

-  'i|  empirical  relation  (12:34): 


t 


d 


b 

I  dt 


(5.1) 
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where  1  is  the  power  per  unit  area,  ^t  is  the  'iuration  of 
the  write  pulse,  and  b  is  a  proportionality  constant.  this 
the  memory  tiae  of  the  device  is  inversely  proportional  to 
the  incident  energy. 

In  the  aajority  of  this  study  it  was  asauaed  that  the 
injection  current  across  the  illuminated  face  was  zero.  In 
selective  tests,  the  boundary  condition  was  aodified  to 
allow  all  of  the  charge  escaping  through  the  dark  face  to 
re-enter  through  the  illuminated  face.  This  is  equivalent 
to  assuming  the  barrier  potential  equal  to  zero.  When  this 
was  done  in  the  aiaple  case  of  the  spot,  the  magnitude  of 
the  integrated  transverse  electric  field  was  seen  to  drop  by 
at  least  three  orders  of  magnitude.  This  suggests  that  the 
injection  current  is  a-major  contributor  to  the  optical 
erasure . 

If  we  assume  a  barrier  potential  at  the  illuminated 
face,  we  would  expect  a  large  magnitude  electric 
field  at  the  illuminated  face  to  cause  a  larger  injection 
current  than  a  smaller  magnitude  field.  Thus,  the  larger 
the  value  of  the  electric  field  at  the  crystal,  the  shorter 
the  lifetime  of  the  image.  The  strength  of  the  field  at  the 
illuminated  face  of  the  crystal  is  proportional  to  the  total 
charge  contained  within  the  crystal.  The  rate  of  formation 
of  total  charge  is  proportional  to  the  intensity  of  the 
write  beam.  The  higher  the  intensity  of  the  write  beam,  the 
larger  the  magnitude  of  the  electric  field  at  the 
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i  1  i  urn  i  mt  ».#»<!  face  of  the  crystal,  and  the  shorter  the  decay 
Mse  of  the  image,  t^. 

This  suggests  two  areas  of  further  research.  One  is 
the  aore  accurate  aodeling  of  the  injection  current.  ThLS 
will  involve  research  into  the  nature  of  the  barrier 
potential  at  the  illuminated  face  of  the  crystal.  The 
second  area  of  research  is  in  the  area  of  increasing  the 
memory  t.iae  of  the  PRIZ.  Both  theoretical  and  experiaental 
work  should  be  pursued  with  a  " semi -conduct ing"  PRIZ,  where 
the  charge  is  allowed  to  escape  freely  from  the  dark  face  of 
the  crystal,  and  a  dielectric  is  used  to  prevent  electrons 
from  entering  the  illuminated  face  of  the  crystal.  It  is 
postulated  that  significant  memory  times  could  be  achieved 
even  with  fairly  strong  write  beams  in  this  manner. 

The  only  solution  to  the  actual  integrated  transverse 
field  presented  in  this  study  is  the  case  of  a  simple  spot. 
In  addition,  the  transverse  field  internal  to  the  write  beam 
is  neglected.  More  work  on  the  solution  to  the  transverse 
field  in  more  complex  cases  is  needed,  in  addition  to 
modeling  the  fields  effects  on  the  electrons  in  the 
conduction  band.  A  good  starting  point,  using  the  charge 
distribution  generated  in  this  study,  would  be  to  modify  the 
results  of  Owechko  and  Tanguay  (ref.  14)  to  model  much  more 
complex  images. 
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